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SUMMARY 


A  series  of  metal  ana  dielectric  targets  exhibiting  fractal  characteristics  are  produced 
having  a  range  of  surface  roughnesses  between  about  1pm  and  16pm.  They  are  all 
sufficiently  rough  that  only  diffuse  reflection  was  obi* 'fed  from  them  when  illuminated 
with  a  He-Ne  laser  beam.  The  intensity  scattered  in  "he  specular  direction  has  been 
measured  as  a  function  of  the  surface  roughness  for  near  normal  (ftj  =  20  )  and  near 
grazing  ( Pj  =  85  )  incidence.  The  specular  direction  intensity  at  first  falls  rapidly  with 
increasing  surface  roughne  s,  but  then  saturates  at  high  values  of  the  surface  roughness 
Good  agreement  is  obtained  between  the  experimentally  observed  variation  with  surface 
roughness  and  that  predicted  theoretically  using  a  fractal  surface  scattering  model  and 
surface  parameters  measured  using  a  profile  instrument.  A  simple  graphical  model  is  used 
to  demonstrate  qualitatively  .lie  trends  observed. 
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1  Introduction 

The  aim  of  this  Memorandum  is  to  present  experimental  results  of  He-Ne  laser  scattering 
near  grazing  incidence  from  surfaces  that  are  very  rough  compared  to  the  laser 
wavelength;  for  comparison  results  are  also  presented  for  near  normal  incidence.  Only 
forward  scattering  is  considered.  Future  work  will  concentrate  on  grazing  angle  backscatter 
geometries  because  of  their  relevance  to  high  resolution  radar  scattering  from  the  sea 
surface. 

Little  data  exists  on  either  forward  grazing  angle  scattering  or  forward  scattering  from  very- 
rough  surfaces  Work  in  the  optical  field  has  tended  to  concentrate  on  surfaces  having 
roughness  ranging  from  considerably  less  than,  to  slightly  greater  than  the  optical  probing 
wavelength.  This  regime  is  of  limited  usefulness  however  for  trying  to  elucidate  some  of 
the  mechanisms  responsible  for  observations  made  of  radar  scattering  from  the  sea-surface 
Examples  of  these  are  the  ratio  of  vertical  to  horizontal  polarised  returns  in  backscatter 
and  the  variation  of  the  returns  with  sea-surface  roughness  in  both  forward  and 

backscatter  geometries. 

There  is  at  present  no  universally  accepted  sea-surface  model.  However,  one  that  has 
had  some  success  in  allowing  some  aspects  of  radar  scattering  from  the  sea  to  be 
predicted  is  the  composite  model'' );  this  assumes  the  sea  is  composed  of  two  different 

scales  of  waves  which  affect  the  radar  cross-section  in  different  ways.  A  spectrum  of 
small  capillary  waves  with  amplitudes  less  than  the  radar  wavelength  is  assumed  which 
cause  Bragg  resonant  scattering.  These  small  capillary  waves  are  in  turn  assumed  to  be 

distributed  on  larger  facets  which  are  tilted  from  the  horizontal  to  represent  the  larger 

gravity  waves  or  swell  Unfortunately  this  model  fails  to  predict  the  ratio  of  vertical  to 
horizontally  polarised  radar  cross-sections  observed  from  backscatter  measurements'2). 
Recently  a  sea-surface  model  based  on  a  fractal  geometry  has  been  proposed  This  is 
discussed  in  more  detail  in  references  3  and  4.  At  present  it  appears  that  a  fractal  index 
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i  -v  2,  corresponding  to  a  marginal  fractal  may  be  realistic,  at  least  for  some  sea 
surfaces. 

In  this  Memorandum  a  series  of  surfaces  having  fractal  properties  with  measured  surface 
roughnesses  ranging  from  just  greater  than  the  optical  wavelength  (X  =  0.6328  pm)  to 
approximately  25X  are  investigated.  The  surfaces  are  described  in  detail  in  section  2  and 
the  scattering  measurements  in  sections  3,  4  and  5.  A  comparison  of  the  measured  mean 
scattered  intensity  with  theoretical  predictions  is  made.  The  results  of  some  simple 
graphical  simulation  are  presented  in  section  6  as  an  aid  to  understanding  the  observed 
trends.  Future  work  will  attempt  to  produce  two  scale  surfaces  as  well  as  fractal  ones 
and  determine  their  different  scattering  characteristics  in  both  backscatter  and  forward 
scatter  geometries. 

2.  Surface  Preparation 

A  series  of  stainless  steel  targets  of  varying  roughnesses  were  prepared  by  randomly 
pressing  an  engraving  tool  onto  the  metal.  This  technique  has  been  previously  shownl^) 
to  produce  a  good  approximation  to  a  fractal  surface  over  a  wide  range  of  scale  sizes 
when  applied  to  Germanium.  By  varying  the  pressure  on  the  engraving  tool  the  surface 
roughness  could  be  varied.  Unfortunately  this  technique  naturally  resulted  in  the  surface 
height  autocorrelation  function  p  increasing  with  the  surface  roughness.  No  simple  method 
could  be  found  of  independently  varying  the  surface  roughness  and  p  over  the  range  of 
roughnesses  required. 

Surface  height  profiles  along  one-dimensional  scans  across  the  surfaces  were  recorded  using 
a  'Talysurf'  stylus  instrument.  From  these  the  surface  stiuctuu.  function,  defined  as 

S(6)  =  <[h(x)  -  h(x  +  6}2>  (1) 

where  h  is  the  local  surface  height  at  position  x,  was  determined.  For  a  fractal  surface 
this  is  related  to  the  two  fractal  parameters  i  and  L  that  completely  describe  it  by 

S(i)  =  L2->-  i  8  i"  (2) 

The  fractal  index  v  is  related  to  the  Hausdorff-Besicovitch  dimension  D  of  a  section 
profile  of  the  surface  by  >•  =  2(2-D)  and  L  is  the  topothesey;  this  is  the  length  over 
which  the  rms  slope  difference  is  unity.  From  (2)  it  follows  that  a  plot  of  log  S(6) 
verses  log  8  allows  the  determination  of  both  »•  and  L.  A  typical  structure  function  plot 
is  shown  in  Figure  1 ;  the  straight  line  has  been  drawn  through  the  experimental  points  for 
reference.  It  is  obvious  that  the  plot  is  linear  over  a  range  of  scale  sizes  from  3pm, 
which  is  the  lower  limit  of  resolution  of  the  'Talysurf  up  to  ~  70pm.  The  limited 
horizontal  resolution  of  the  stylus  instrument  unfortunately  precludes  characterising  the 
surface  at  scale  sizes  comparable  to  the  optical  probing  wavelength  of  0.633  pm. 
However  there  is  no  a-priori  reason  for  assuming  that  the  fractal  regime  indicated  by  the 
linear  form  of  the  log-log  structure  function  does  not  extend  down  to  this  region;  the 
scattering  data  presented  later  re-enforces  this  view.  The  slope  of  the  straight  line  region 
gives  i'  and  the  intercept  allows  L  to  be  determined. 

A  surface  roughness  parameter  is  invariably  used  when  describing  rough  surfaces 
Unfortunately  such  parameters  are  obviously  beset  with  difficulties  when  applied  to  fractal 
surfaces.  For  example,  the  measured  surface  roughness  is  a  function  of  the  length  over 
which  it  is  measured.  Nevertheless,  because  of  the  need  to  make  comparisons  with  sea 
scattering  data  etc  where  the  surface  is  described  in  such  terms,  the  arithmetic  average 
roughness  Ra  was  measured  for  each  surface  with  a  small  portable  stylus  instrument  (Rank, 
Taylor,  Hobwin  Surtronic);  this  measured  Ra  over  a  5  mm  length  of  surface  and  used  a 
0.8  mm  cut-off.  A  relation  between  Ra  and  the  fractal  parameters  r  and  L  is  described 
in  Section  5. 
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1  he  arithmetic  average  lcughness  is  defined  for  a  profile  length  L  (5  mm  in  these 
measurements)  as 


R 


a 


L 

izi  dx 


where  z  is  measured  from  the  mean  level.  For  a  symmetrical  Gaussian  distribution  of 
surface  heights  Ra  is  related  to  the  rms  surface  roughness  a  by 


R 
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Surfaces  were  produced  with  Ra  values  of  0.92pm,  2.3pm,  4.7pm,  6.3pm,  6.6pm,  10.2pm 
and  IS. 5pm,  covering  the  range  1 .4  ^  Ra/X  <  24.5.  For  all  the  samplts  the  fractal  index 
r  was  approximately  constant  and  equal  to  1.66  t  0.04.  The  topothesey  however  varied 
over  quite  a  wide  range,  from  2  x  10-^  pm  to  6.3  x  10~^  pm.  All  the  samples 
possessed  an  isotropic  surface. 

The  measured  surface  height  histograms  for  surfaces  with  Ra  =  2.3pm,  10.2pm  and 
15.5pm  are  shown  in  Figures  2(a),  2(b)  and  2(c)  respectively.  The  vertical  axis  is  the 
number  of  points  having  heights  (in  arbitrary  units)  in  the  range  given  by  the  horizontal 
axis.  Also  shown  for  each  one  is  a  Gaussian  distribution  calculated  using  the  measured 
mean  and  standard  deviation.  From  these,  which  are  typical  of  all  the  surface  height 
histograms,  it  can  be  seen  that  at  best  the  histograms  are  only  very  approximately  fitted 
by  Gaussian  distributions.  The  measured  Skewness^)  and  KurtosisW  js  shown  for  each 
one  in  the  figures.  For  comparison  an  attempt  to  fit  the  histograms  to  a 
Gamma-distribution  was  also  made.  For  this  the  height  distribution  is  given  by 


...  r  m  tm  ,m-l  exp  I-  mh/<h>] 

p(h)  -  l  1  h  - roo - 

where  the  fitting  parameter  m  was  determined  from  the  measured  values  of  <h^>  and 
<h>^  as 


<h2> 


<h> 


+ 


1 
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As  can  be  seen  from  Figures  2(a),  2(b)  and  2(c)  the  resulting  fit  to  the  measured  data 
was  if  anything,  slightly  worse  than  that  obtained  using  a  Gaussian  distribution. 

For  comparison  purposes  a  few  surfaces  were  prepared  using  dielectric  (perspex;  rather 
than  a  metal.  Three  surfaces  having  Ra  values  of  0.75  pm,  5.8pm  and  14.2  pm  were 
made  The  fractal  index  r  was  rather  higher  in  this  case,  being  approximately  1  84  (t 
0.03).  L  ranged  from  1.3  x  10-*  to  4.2  x  10_!^  pm. 

3.  Experiment 

The  experimental  arrangement  is  shown  schematically  in  Figure  3.  The  output  from  a 
polarised  He-Ne  laser  was  passed  through  a  chopper,  scattered  from  the  target  and 
detected  by  a  3  mm  wide  by  9  mm  long  photodiode.  The  output  from  the  detector, 
togther  with  the  chopper  reference  signal  was  connected  to  a  Brookdeal  5209  Lock-in 
amplifier.  All  measured  intensities  are  quoted  in  micro-volts.  The  detector  was  mounted 
on  an  arm  that  replaced  the  telescope  on  a  divided  circle  spectrometer,  the  collimator  arm 
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having  been  removed.  The  target  was  mounted  on  the  prism  table  at  the  centre  of  the 
spectrometer.  The  detector  to  centre  of  target  distance  was  ~  370  mm 

The  laser  beam  diameter  was  3.5  mm  at  the  target  position.  Angles  of  incidence  Pj  of 
greater  than  -  87’  led  to  illuminated  lengths  comparable  to  the  target  size  of  40  mm  and 
so  was  limited  to  a  maximum  of  85’.  Because  of  the  relatively  small  illuminated  size 
near  normal  incidence  the  experiments  were  repeated  several  times  with  the  target  moved 
between  each  measurement  and  the  results  averaged  to  ensure  that  a  large  proportion  of 
the  target  was  sampled.  The  targets  were  also  rotated  to  ensure  they  were  isotropic. 
Speckle  averaging  was  automatically  achieved  by  the  use  of  a  relatively  large  detector 

aperture  (3  mm  x  9  mm),  the  typical  speckle  size  being  of  order  70  pm. 

Measurements  were  carried  out  with  the  incident  radiation  either  p  or  s  polarised. 
P-polarised  light  has  its  electric  field  vector  parallel  to  the  plane  of  incidence  and 

corresponds  to  V-polarised  radiation  in  radar  terminology;  s-polarised  light  (H-polarised  in 
radar)  has  its  electric  field  vector  perpendicular  to  the  plane  of  incidence.  Changes  fiom 
one  state  of  polarisation  to  the  other  were  accomplished  by  rotating  the  cylindrically 
symmetrical  laser  in  its  mounting  The  detection  system  was  made  polarisation  sensitive 
by  mounting  a  .v/iatable  sheet  polariser/analyser  of  extinction  ratio  approximately  400.1  in 
front  of  the  photodiode. 

4.  Diffuse  Reflection 

To  verify  that  the  metal  targeb  were  actually  acting  as  fractal  scatterers  the  mean 
scattered  intensity  as  a  function  of  angle  <1(P)>  was  measured  for  one  particular  target 
for  angles  of  incidence  Pj  of  7.0  and  85  and  compared  with  that  predicted  from  theory 
using  'Tallysurf  measured  values  of  v  and  L.  The  target  with  Ra  =  2.3  pm  was  used; 

for  this  surface  r  -  1.63  and  L  =  5  x  10~4pm. 

For  a  perfectly  conducting  isotropic  two-dimensional  fractal  surface  the  diffusely  scattered 
component  is  given  by(6) 

2  r  a  ■>  2-i  <■ 

<I(0)>  <*  F  (0)  j  r  J  ‘ krC ( 0 ) j  exp  [  - ;  k  r*(0)  L  r  j  dr  (3, 

o 


where  various  inessential  angle-independent  terms  h3ve  been  omitted  F,  G  and  f  are 
angle  dependent  terms  given  by 


1  +  cos  0  j  cos  (>2  -  sin  P  j  sin  P^  cos 

Ftp)  -  - - - - - 

cos  Pj  +  cos  Pj 

2  2  * 

G ( 0 )  -  [(sin  Sj  -  sin  0^  cos  0^)  +  sin  0^  sin 

f (0)  -  cos  0j  +  cos  0^ 


(4) 


(5) 

(6) 


where  0j ,  Pj,  and  Pj  are  as  defined  in  Figure  4  For  all  the  measurements  described  in 
this  Nitiiioidiiuuni  03  =  0 

Equation  (3)  should  strictly  be  limited  to  situations  where  the  Kirchoff  approximation  is 

valid  i.e  to  situations  in  which  small  angle  scalar  diffraction  theory  is  valid.  This  of 

course  is  questionable  for  fractals  because  they  include  scale  sizes  less  than  the  wavelength 
of  the  incident  radiation.  Provided  however  the  scattering  power  of  the  inner  scale  sizes 
is  small,  by  which  is  meant  their  rms  'height'  is  much  smaller  than  a  wavelength,  then 

the  basic  scattering  characteristics  of  the  fractal  model  will  be  unaffected  in  near  specular 


4 


directions.  Small  scale  roughness  will  however  contribute  to  large  angle  scattering  and  so 
some  doubt  exists  as  to  the  validity  of  Fraunhofer  large  angle  scattering  results  obtained 
using  a  fractal  model.  There  are  also  a  number  of  other  approximations  implicit  in  using 
(3)  at  large  scattering  angles.  These  include  neglect  of  the  variation  of  reflectivity  with 
angle,  shadowing  and  polarisation  effects  which  are  obviously  not  included  in  the  scalar 

treatment  used  here.  However,  in  spite  of  these  omissions  there  is  an  accumulating  body 
of  experimental  evidence  that  exhibits  trends  consistent  with  those  predicted  by  (3)  even 
from  multi-scale  surfaces  at  relatively  large  scattering  angles^'^.  Near  grazing  angle 
results  have  not  however  previously  been  published. 

The  variation  of  <I(tf)>  with  scattering  angle  8  for  =  20'  and  85"  is  shown  in 

Figures  5  and  6  respectively;  angles  are  measured  from  the  normal  to  the  mean  target 

surface.  Also  shown  are  the  curves  calculated  from  (3)  using  the  measured  values  of  v 

ana  L.  In  both  figures  the  theoretical  curve  has  been  normalised  to  the  mean  of  the 

measured  p  and  s-polarised  returns  at  6  =  38'.  The  good  agreement  in  shape  between 
the  measured  and  calculated  curves  is  readily  apparent.  In  both  cases  the  calculated  curve 
underestimates  the  measured  .«lues  near  the  specular  direction,  and  this  is  discussed  more 
fully  later.  It  also  fails  at  angles  above  about  85'  for  the  8t  =  20’  case  due  to  the  lack 
of  any  shadowing  correction  in  (3).  A  simple  shadowing  correction  that  has  been 

successfully  used  before  is  outlined  in  reference  6.  The  good  agreement  between  the  form 
of  the  calculated  and  measured  angular  variations  of  the  scattered  diffuse  intensity,  even  at 
near  grazing  incidence  lends  support  to  the  contention  that  the  metal  surfaces  do  indeed 

behave  as  fractal  scatterers  with  the  measured  values  of  i-  and  L. 

5.  Specular  Reflection 

The  variation  of  the  mean  intensity  scattered  in  the  specular  direction  <I(0S)>  with 

surface  roughness  Ra  is  shown  in  Figure  7  for  8 j  =  20’  and  Figure  8  for  =  85'. 

There  is  very  little  difference  between  the  p  and  s-polarised  results  for  =  20  ;  nearer 
grazing  incidence  (6;  =  85  )  however  the  s-polarised  signal  is  consistently  higher  than  the 
p-polarised  one  for  all  values  of  Ra  by  a  factor  of  3  to  4.  This  is  purely  a 

manifestation  of  the  different  behaviour  of  the  Fresnel  reflection  coefficients  for  p  and 
s-polarised  radiation  with  angle,  and  is  determined  not  by  the  surface  roughness  but  by 
the  complex  refractive  index  of  the  metal  target.  This  is  illustrated  in  Figure  9  which 

shows  the  variation  of  the  specular  signal  with  angle  of  incidence  for  the  smoothest  metal 
surface  investigated  (Ra  =  0.9  pm).  It  shows  there  is  very  little  difference  between  the  p 
and  s-polarised  signals  at  8 j  =  20',  but  at  85  the  s-polarised  signal  is  about  3.5  times 

that  of  the  p-polarised  one.  As  grazing  incidence  is  approached  the  p  and  s-polarised 
returns  became  closer  to  each  other 


The  results  shown  in  Figures  7  and  8  for  both  angles  of  incidence  show  that  the  reduction 
of  specular  direction  intensity  with  increasing  surface  roughness  levels-off  at  large  values  of 
1*3,  the  effect  being  more  noticeable  for  the  large  angle  of  incidence  case.  Some  simple 
graphical  modelling  to  show  the  physical  basis  for  this  behaviour  is  presented  in  Section  6. 

For  all  the  targets  investigated  in  this  work  the  surface  constitutes  a  very  rough  target 
((4t/X)<t  cos  8 j  >  1)  even  at  angles  of  incidence  of  85'.  Consequently  a  true  specular 
component  is  unlikely;  the  reflected  radiation  in  the  specular  direction  should  be  diffusely 
reflected  Therefore  equation  (3)  may  be  expected  to  be  appropriate  and  can  be 
evaluated  analytically  in  the  specular  direction  for  arbitrary  values  of  v  and  L  as 


< I ( P  )>  a 

s 


F2(ft) 


— n — 

.■[Jk  f  (f) 


2-r  2/r 


(7) 


The  values  calculated  using  this  expression  are  shown  superimposed  on  the  measured  values 
in  figures  7  and  8  and  marked  'Theory'.  As  with  the  calculations  of  the  angular 
variation  of  scattered  radiation  shown  in  Figures  5  and  6,  the  results  have  been  scaled  by 
an  appropriate  factor  to  bring  them  into  agreement  with  the  measured  values  This  is 
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necessary  because  of  the  omission  of  various  inessential  angle  independent  factors  in  (3). 
The  reason  for  the  difference  in  scaling  factor  between  the  specular  and  non-specular 
cases  remains  obscure.  It  is  obvious  from  Figs.  5  and  6  that  the  scaling  factor  required  to 
fit  the  data  away  from  the  specular  direction  underestimates  the  specular  direction  values. 
However,  the  excellent  agreement  shown  in  Figs  7  and  8  between  the  measured  specular 

direction  values  and  those  calculated  using  an  appropriate  scaling  factor  strongly  suggests 

that  (7)  is  correct.  Work  by  Kim  et  al®  indicates  that  for  diffusers  with  Gamma 

distributed  surface  profiles  the  specular  component  is  much  higher  than  expected 
Experimentally  they  observed  a  specular  component  with  He-N’e  laser  radiation  when  the 
standard  deviation  of  the  surface  height  fluctuations  was  as  high  as  2pm  from  a  target 

possessing  a  negative  exponential  surface  height  distribution.  However  the  surface  height 

histograms  shown  in  Figure  2  do  not  appear  to  be  obviously  Gamma -distributed  and  so 
the  reason  for  this  behaviour  remains  unresolved  at  present. 

Traditionally  measurements  of  surface  structure  are  made  in  terms  of  the  surface  roughness 
a  and  the  surface  height  correlation  length  C,  not  in  terms  of  fractal  parameters. 

Because  a  simple  finite  scale  model  behaves  like  a  fractal  at  small  displacements  a 

correspondence  may  be  drawn  between  such  models  and  fractal  ones  For  any  model 
possessing  a  Eorentzian  spectrum  and  surface  rough  enough  to  produce  large  phase 

fluctuations  on  reflection,  it  can  be  shown  that  for  two-dimensional  isotropic  surfaces 


Assuming  a  Gaussian  height  distribution  a  can  be  expressed  in  terms  of  Ra  as 
c  -  125  Ra.  For  the  surfaces  studied  in  this  work  the  surface  height  correlation  length 
f  varies  linearly  with  Ra  as  shown  in  Figure  10.  i.e.  C  =  4Ra  +  B.  Substituting  this  in 
(8)  leads  to  the  following  expression  for  the  topothesey  in  terms  of  r  and  Ra  as 


1  2 
„  a  ( ^  "  >  )  1 1  i  <  %  2  - 


From  Fig  10,  A  =  5.5  and  B  =  12  pm.  The  variation  of  L  with  Ra  calculated  from  (9) 
is  shown  in  Figure  11  and  is  in  reasonable  agreement  with  the  values  determined  from  the 
structure  function  plots. 

A  similar  variation  of  the  specular  intensity  with  surface  roughness  * o  that  found  from  the 
metal  targets  was  observed  using  the  perspex  ones;  these  are  shown  in  Figure  12.  For 

ease  of  following  the  trends  a  continuous  curve  has  been  drawn  through  each  set  of 
results.  As  expected  the  biggest  difference  between  the  results  obtained  using  the  metal 

and  dielectric  targets  is  in  terms  of  the  degree  of  cross  depolarisation  observed.  This  is 
defined  (as  a  percentage)  as  the  ratio  of  the  signal  delected  with  the  polariser/analyser 

orthogonal  to  that  emitted  by  the  laser,  to  that  with  it  parallel  to  the  emitted  polarisation 

direction  Some  results  for  the  metal  target  are  shown  in  Figure  13;  to  make  the  plots 
easier  to  read  the  measurements  at  discrete  angles  have  been  replaced  by  continuous 
curves.  Depolarisation  is  shown  plotted  against  angle  for  two  angles  of  incidence  (20 ’  and 
85  )  for  s  and  p-polarisation  of  the  laser.  As  seen,  the  degree  of  depolarisation  is  quite 
small,  ranging  from  1%  or  less  near  the  specular  direction  to  a  maximum  of  about  10% 
well  away  from  it.  The  values  obtained  well  away  from  the  specular  direction  are 
possibly  an  overestimate  as  the  depolarised  component  was  very  small  and  comparable  to 
the  system  noise  level  For  the  perspex  targets  the  degree  of  depolarisation  near  the 
specular  direction  was  again  of  order  1%  or  less,  but  in  this  case  rose  to  very  high  values 
away  from  it  as  shown  in  Figure  14.  For  example  for  p-polarised  light  the  depolarisation 
approached  100%  as  expected  at  the  Brewster  angle  ("-  56  ), 


Work  in  the  microwave  radar  field  on  near  grazing  angle  forward  scatter  had  led  to  the 
suggestion'®  that  a  form  of  roughness  parameter  defined  as  (<j  cos  0j/X)  may  be  useful 
The  suggestion  is  that  if  the  surface  roughness  rr  and  the  angle  of  incidence  S,  are  varied 
so  as  to  maintain  a  cos  dj  constant,  the  scattering  signal  should  remains  constant  unless 
shadowing  is  important.  It  may  therefore  be  used  to  estimate  the  importance  of 
shadowing  A  brief  attempt  has  been  made  to  check  this  idea.  Figure  15  is  a  plot  of 

the  specular  intensity  from  the  metal  targets  as  a  function  of  Rj  for  Ra  cos  0j/x  =  1 .45 
The  lowest  Rj  surface  corresponds  to  an  angle  of  incidence  of  20’  and  the  highest  one  to 
dj  =  82".  The  same  data  is  shown  plotted  as  a  function  of  the  angle  of  incidence  in 
Figure  16  Obviously  the  specular  direction  signal  is  far  from  constant  as  the  angle  of 
incidence  and  surface  roughness  is  varied.  Although  the  data  is  very  limited  and  any 
conclusions  drawn  must  be  tentative  ones.  Figure  16  suggests  that  the  specular  direction 
signal  may  be  relatively  constant  from  low  angles  of  incidence,  up  until  perhaps  60  -70  , 
beyond  this  the  signal  fails  off  markedly.  Similar  results  were  obtained  for  Racos0/X  = 
2  1.  Care  must  be  exercised  in  transferring  these  tentative  findings  directly  to  the  radar 
situation  because  of  differences  in  radar  transmitter  and  receiver  antennae  gains,  and 
optical  detector  field  of  view  compared  to  the  illuminated  target  region.  For  example,  as 
grazing  angle  is  approached  the  illuminated  patch  becomes  extremely  long  and  may  well 
exceed  the  receiver  field  of  view  More  work  is  planned  in  this  area  to  clarify  some  of 

the  uncertainties. 

6.  Graphical  Simulation 

As  a  simple  aid  to  understanding  the  reasons  for  the  fall-off  and  saturation  of  the 
specular  direction  intensity  with  increasing  surface  roughness  a  simple  graphical  approach 
was  adopted.  A  rough  surface  profile  was  generated  on  a  computer  from  a  sum  of 
sinusoidal  waves  and  the  effect  of  varying  the  roughness  monitored.  Initially  a  sum  of 
equally  phased  large  scale-size  waves  of  the  same  amplitude  was  used,  figure  17  shows 
such  a  profile.  The  amplitude  A  is  given  by 


A  =  150  +  o  (sin  3x  +  sin  1.5x  +  sinx  *  sin  0  75x  +  sin  0.5x)  (10) 

A  histogram  of  the  surface  height  distribution  is  shown  in  figure  18  for  o  =  10  Also 

shown  is  a  Gaussian  distribution  of  the  same  mean  and  standard  deviation  as  the  measured 
one.  Because  of  the  equality  of  the  phases  of  the  component  waves  the  surface  heights 
are  not  truly  Gaussian  distributed,  this  is  also  true  for  the  actual  metal  and  perspex 

surfaces  described  earlier  (Figure  2).  The  e"'  correlation  length  was  46.  As  expected, 

the  structure  function  plot  for  this  surface  indicated  a  fractal  index  of  ~  2  and  a 

topothesey  of  approximately  zero  as  shown  in  Figure  19. 

As  specular  reflection  only  occurs  from  small  regions  or  facets  having  zero  slope,  the 

specular  intensity  was  taken  as  the  number  of  points  on  the  profile  having  a  slope  less 

than  some  arbitrary  value;  in  this  case  less  than  one  degree  Using  this  criterion  the 
effect  of  varying  o  and  hence  the  surface  roughness  on  tne  specular  intensity  is  shown  in 
Figure  20  for  three  angles  of  incidence  It  is  immediately  obvious  that  tht  same  genera! 
trend  of  an  initially  rapid  fall-off  of  intensity  with  increasing  roughness,  as  was  seen  with 
the  actual  rough  surface  measurements  is  observed  The  reduction  in  specular  intensity 
with  increasing  a  is  greatest  for  the  largest  angles  of  incidence,  and  is  due  to  shadowing 
The  effect  on  the  surface  slope  distribution  of  increasing  the  surface  roughness  is  shown  in 
Figure  21,  and  demonstrates  the  very  large  reduction  in  the  number  of  low  slopes  as  the 
surface  roughness  is  increased 

Another  way  of  dramatically  affecting  the  surface  slope  distribution  is  to  add  a  high 

frequency  component  to  the  wave  structure.  This  is  also  a  possible  moJei  for  waves  on 

the  sea  surface.  To  this  end  the  wave  given  by  (10)  was  modified  by  the  addition  of 
components  of  ten  times  the  spatial  frequency  of  those  in  (10)  to 
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A(x)  =  150  +  o  (sin  3x  +  sin  1.5x  +  sin  x  +  sin  0  ?5x  +  sin  0  5x) 

+  i«/4  (sin  30x  +  sin  l5x  +  sin  lOx  +  sin  7  5x  *  sin  5x)  (11) 

Such  a  surface  with  n  =  10  is  shown  in  Figure  22  The  much  more  irregular  appearance 

of  this  surface  profile  compared  to  the  original  shown  in  Figure  17  is  immediately 
apparent.  The  effect  of  adding  the  high  frequency  components  on  the  structure  function 

of  the  original  low  frequency  waves  is  shown  in  Figure  19.  The  normalised  surface  height 
autocorrelation  function  of  the  large  scale  waves,  small  scale  ones  and  their  sum  is  shown 
in  figure  23  The  effect  on  the  surface  slope  distribution  is  given  in  Figure  24.  and 

shows  a  similar  effect  to  increasing  the  roughness  of  the  large  scale  only  waves  The 
variation  of  specular  intensity  with  surface  roughness  is  shown  in  Figure  25.  At  normal 
incidence  the  effect  of  adding  the  high  frequency  waves  is  to  significantly  reduce  the 
signal  At  larger  angles  of  incidence  the  effect  is  very  dramatic  and  illustrates  quite  well 
the  major  influence  that  small  amplitude  high  frequency  waves  superimposed  on  large  scaie 
waves  have  on  the  specular  direction  scattering  properties. 

7  Conclusions 


The  experimental  measurements  made  utilising  very  rough  surfaces  that  exhibit  fractal 

properties  show  that  the  intensity  scattered  in  the  specular  direction  decreases  with 
increasing  roughness,  and  m  fact  saturates  at  sufficiently  high  values  of  roughness  This  is 
in  qualitative  agreement  with  observations  obtained  using  microwave  forward  scattering  from 

the  sea  surface  close  to  grazing  incidence  Elementary  graphical  analysis  suggests  this  is 

due  to  a  reduction  in  the  probability  of  finding  zero  slope  regions  on  the  profile 
Although  this  concept  is  strictly  only  valid  for  non-fractal  surfaces,  in  practice  it  mav  be 
expected  to  give  a  reasonable  guide  to  the  behaviour  for  values  of  the  fractal  index  of 
order  two.  At  angles  of  incidence  grcaier  than  about  sixty  degrees  shadowing  becomes 
increasingly  important  The  only  signficant  difference  observed  between  relatively  good 

conductors  and  dielectrics  was  in  the  degree  of  cross  depolarisation  produceJ 

Future  experiments  are  planned  in  which  both  the  targets  characterised  in  this  work,  and 

hopefully  targets  containing  two  dominant  scale  sizes  wall  be  investigated  in  a 

backscattering  configuration  This  will  allow  comparisons  to  be  made  with  microwave  ladat 
measurements  of  internal  wave  wakes 
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FIG  2.  SURFACE  HEIGHT  HISTOGRAMS 


FIG  3.  EXPERIMENTAL  ARRANGEMENT 


FIG  4.  GEOMETRY  FOR  SCATTERING  FROM  A  ROUGH  SURFACE 
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FIG  5.  VARIATION  OF  SCATTERED  INTENSITY  WITH  ANGLE  FOR  B,  =  20” 
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FIG  7.  VARIATION  F  SPECULAR  INTENSITY  WITH  SURFACE  ROUGHNESS  FOR 
0j  =  20'.  THEORY  IS  CALCULATED  FROM  (7) 
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FIG  9.  VARIATION  OF  SPECULAR  REFLECTION  WITH  ANGLE  OF  INCIDENCE 
FROM  SMOOTHEST  METAL  SURFACE 
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FIG  1 1 .  MEASURED  AND  CALCULATED  VALUES  OF  THE  TOPOTHESEY  L 
AS  A  FUNCTION  OF  SURFACE  ROUGHNESS 
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FIG  12.  VARIATION  OF  SPECULAR  DIRECTION  INTENSITY  WITH  SURFACE 
ROUGHNESS  USING  PERSPEX  TARGETS.  THE  CONTINUOUS  CURVES 
ARE  PURELY  AN  AID  TO  VIEWING  THE  GRAPH 
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FIG  13  DEPOLARISATION  FROM  METAL  TARGETS  FIGURES  ON  THE 
CURVES  ARE  THE  SURFACE  ROUGHNESS  Ra 
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FIG  15.  VARIATION  OF  SPECULAR  DIRECTION  SIGNAL  WITH 
SURFACE  ROUGHNESS  FOR  (Ra  cos  9;)/?.  =  1.45 
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FIG  21 .  EFFECT  OF  SURFACE  ROUGHNESS  ON  SLOPE  DISTRIBUTION  OF 
LARGE  SCALE  SIZES  ONLY  SURFACE 


FIG  23.  NORMALISED  SURFACE  HEIGHT  SPATIAL  AUTOCORRELATION  FUNCTIONS 
FOR  SURFACES  CONTAINING  LARGE  SCALE  SIZES  ONLY  (10),  SMALL  SCALE  SIZES  ONLY, 
AND  BOTH  LARGE  AND  SMALL  SCALE  SIZES  (11).  FOR  ALL  SURFACES  a  =  10. 
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FIG  24.  EFFECT  ON  SURFACE  SLOPE  DISTRIBUTION  OF  ADDING  HIGH 
SPATIAL  FREQUENCY  COMPONENTS 

a)  LARGE  SCALE  SIZES  ONLY  (a  =  10) 

b)  LARGE  AND  SMALL  SCALE  SIZES  (a  =  10) 
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